Adiabatic dynamics in passage across quantum critical lines and gapless phases 
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It is well known that the dynamics of a quantum system is always non-adiabatic in passage 
through a quantum critical point and the defect density in the final state following a quench shows 
a power-law scaling with the rate of quenching. However, we propose here a possible situation where 
the dynamics of a quantum system in passage across quantum critical regions is adiabatic and the 
defect density decays exponentially. This is achieved by incorporating additional interactions which 
lead to quantum critical behavior and gapless phases but do not participate in the time evolution 
of the system. To illustrate the general argument, we study the defect generation in the quantum 
critical dynamics of a spin- 1/2 anisotropic quantum XY spin chain with three spin interactions and 
a linearly driven staggered magnetic field. 
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Understanding the dynamics of a quantum system 
passing through a quantum critical point [H, [Ij has been 
a very active and fascinating area of research in recent 
years (sl-fl^. The dynamical evolution can be initiated 
either by a sudden change of aparameter in the Hamilto- 
nian, called a sudden quench Q , or by a slow (adiabatic) 
quenching of a parameter 0, Q . The related entanglement 
and fidelity properties are also being looked at ^M,- It is 
well known that when a quantum system, initially pre- 
pared in its ground state, is swept adiabatically across a 
quantum critical point, defects are generated in the final 
state of the system due to the critical slowing down which 
forces the dynamics of the system to be non-adiabatic in 
the vicinity of the quantum critical point. In a linear 
passage through an isolated critical point, when a pa- 
rameter {e.g., the magnetic field) is changed in time t 
as t/r, the defect density (n) in the final state scales 
with the quenching rate l/r following the Kibblc-Zurek 
scaling relation [la, 113 given by n ^ ^-dv/('^'^+^) ^^^g 
adiabatic limit (r — )■ cx))[l, [1, [l2|. Here, d is the spa- 
tial dimension, v and z are the correlation length and 
the dynamical exponents, respectively, associated with 
the quantum critical point The experimental verifi- 
cation of such dynamics is now possible by studying the 
dynamics of atoms trapped in optical lattices |18| . 

In this work, we explore the possibility of adiabatic 
dynamics or exponentially decaying defect density even 
in passage through quantum critical points and gapless 
phases. We show that this occurs only in a special situa- 
tion in which an additional term of the Hamiltonian that 
leads to the quantum critical behavior and gapless phases 
does not participate in the dynamics of the system. We il- 
lustrate the general argument using a quantum spin-1/2 
XY chain with multispin interactions and a staggered 
magnetic field which is exactly solvable via the Jordan- 
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Wigner transformation. We thus provide an example of 
a special situation that contradicts the power-law scaling 
of the defect density for a non-random quantum system 

Let us begin the discussion with a d-dinicnsional quan- 
tum Hamiltonian given by 



H = Y,i:\k)H{k)^{k); 

k 

H{k) = ac(fc)i+ (Aa^ + A(fc)CT+ + A*(fc)a_ 



(1) 



where a, c(k) and A(A:) are model dependent functions, 
a'^{i = x,y,z) are Pauli spin matrices, (t+ = ((T^-|-icr^)/2, 
CT_ = (ct^ — iCT^)/2 and 1 denotes the 2x2 identity ma- 
trix. Here, A defines the time-dependent parameter that 
is to be quenched adiabatically. The column vector ?/>(/?) 
defines a two component fermionic operator. Such an ex- 
actly solvable Hamiltonian (with a = 0) is known to rep- 
resent several one- and two-dimensional intcgrable quan- 
tum spin models as the Ising, the XY spin chains 19 1 
and the extended Kitaev model in two dimensions 20 1 
when the spin operators are transformed into spinless 
fermions via the Jordan- Wigner transformation [l9[. In 
the present case, however, the nature of interaction of the 
spin chains we study, necessitates the consideration of a 
two sublattice structure, and hence the fcrmion operator 
4'{k) = (oj:, h-^) where {h-^) denote the Jordan- Wigner 

Fermions for the mode k describing the spins on even 
(odd) sublatticcs. The excitation energy of the Hamilto- 
nian is given by 

(2) 



ef = ac{k) ± + |A|2. 



The phase diagram of the model in the a — A plane can 
be easily obtained. The presence of the additional term 
ac(fc) plays a non-trivial role in determining the phase 
diagram of the model by making excitation energy zero 
for specific values of the parameter and the wavevector 
k. For example, if A(fc) = for the wavevector fco, we 
obtain a critical line (e^ = 0) given by A = ac(fco)- 
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Our interest is in the defect generation when the pa- 
rameter A is quenched in a hnear fashion as t/r from 
t — > — cx) to t — > -|-oo and the system is swept across 
the critical line. Let us assume that at t — > — oo, the 
system is prepared in its ground state such that 

ata^jlj;) — 1 and = for any k. For an adi- 

abatic dynamics, the expected final state is |2^) defined 
as 6t5^|2^) — 1. In course of dynamics, a general state 
describing the system at an instant t can be put in the 
form ^ C^^(t)|lj:) + C.^^{t)\2j^) where the time 

dependent coefficients satisfy the Schrodinger equation 



d 



( m 

V A(fc)* 



l.k 
^2,k 



A{k) \ 

-m J 



(3) 



with h set equal to unity. Also, the time dependent pa- 
rameter X{t) varies as Xot/r and the initial conditions are 
\C^j:it -oo)|2 = 1 and \C^j:it ^ -oo)p = 0. It is 

worth noting that although the term ac{k)l plays a cru- 
cial role in determining the critical line or the gapless re- 
gions, it does not show up in the time-dependent Hamil- 
tonian Hj:{t) which dictates the temporal evolution. This 
is because the identity operator commutes with all the 
other terms of the Hamiltonian at every instant of time. 
The term ac{k)l influences the dynamics only up to a 
phase factor and is hence truly irrelevant in deciding the 
non-adiabatic behavior. In this sense, the time depen- 
dent Hamiltonian Hj:(t) does not capture the passage 
through quantum critical lines and gapless phases gener- 
ated by ac(fc)! The Schrodinger equations (3) describing 
the dynamics of the system effectively boils down to the 
standard Landau-Zener problem (LZ) [2l| of two time- 
dependent levels ±\/A(<P"T1A|^ (not the levels given in 
Eq. 2) approaching each other in a linear fashion with a 
minimum gap 2|A| at time t = 0. The probability of exci- 
tation in the final state is given by Landau Zencr formula 
MM Pfc = |Ci_g(t^oo)|2 =exp(-(^A(fc)2r)/Ao) . 

Let us assume that the parameter A(k) vanishes at a 
quantum critical point for a = as A{k) ^ |fc — fcol^^ 
where fco is the critical wave vector. Noting that in the 
adiabatic limit of large r, only modes close to fcg con- 
tribute, the defect density in the final state is given by 
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{^Trr JBZ 
1 



(27r)'' 
1 

-j-d/2z2 



Pk d'^k 

d'^k exp(-(7r|fc- fcoP^V)/Ao) 



BZ 



(4) 



The scaling of the density of defects hence depends only 
on the exponent 22 as observed previously in references 
PH in the context of quenching through a multicritical 
point. The situation where the parameter A(fc) does 



never vanish during the temporal evolution of H^(t) is 
far more interesting. If the parameter A attains the 
minimum value Aq for some wave vector fco satisfying 
the scaling form A^ = Aq + 5\k — fcoP^^j the nona- 
diabatic transition probability will show an exponen- 
tial behavior and the defect density will scale as n ^ 
exp(— 7rAQT/Ao)/'r'^/^^2^ As discussed already, the scal- 
ing of the defect density satisfy the same scaling form 
for all values of a. This exponential decay of defect even 
in passage through a quantum critical point is the key 
feature of this communication. 

Question is whether it is possible to find a Hamiltonian 
which gets mapped to the Eq. (1). To show this we con- 
sider a spin-1/2 quantum XY spin with a two sublattice 
structure in the presence of a three-spin interaction and 
a staggered magnetic field h given by the Hamiltonian 

i i 

- J2 5Z«2<+1,1 + <2'^f+l,l) - "^3 5I«1<2<+1,1 

i i 

+ <i<2'^f+l,l) - ^3 l]('^^^2'7f+l,l'7f+l,2 

i 

+ '^LCTl+1,1'^1^1,2)' 



(5) 



where i is the site index and the additional subscript 
1(2) defines the odd (even) sublattice. The parameter 
Jl describes the XY interaction between spins on sub- 
lattice 1 and 2 while J2 describes the XY interaction 
between spins on sublattice 2 and 1 such that Ji is not 
necessarily equal to J2. The parameter J3, chosen to be 
positive throughout, denotes the three spin interaction. 
Some variants of the Hamiltonian ([5]) were studied pre- 
viously [2^. This spin chain is exactly solvable in terms 
of Jordan- Wigner fermions [2^ defined on even and odd 



sublattices as a^-^ 



nj<*(-<^l,i)(-<2)J al> and 
b], where af = 2aja, 
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n,<.(-^i,i)(-^i,2)(-<i) 

and (T.fj ^ '2'b\bi — 1. The Fcrmion operators a; and bi 
can be shown to satisfy Fcrmionic anticommutation re- 
lations. 

In terms of the Jordan- Wigner Fermions, the Hamil- 
tonian (5) can be recast in the Fourier space to the form 
given in Hamiltonian (1) with 



a cos fcl 



A -(l+7e-»'^-) 
-(1+76+''^) -A 



(6) 



where ipl, ~ (aj., b\) and we have set A = /i/Ji, a = J3/J1 
and 7 — J2I J\- The excitation energy is obtained as 



I cos fc ± - V'^^ + 1 + 7^ + 27 cos fc 



(7) 



Comparing with the Hamiltonian (1), we also identify 
c(fc) — cosfc and |A(fc)p = 1 -|- 7^ + 27Cosfc. The phase 
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FIG. 1: Phase diagram of the Hamiltonian (5) in the a — A 
plane for 7=1 and 7 = 0.5. Different phases are discussed 
in the main text. The vertical line shows the direction of 
quenching. 



FIG. 2: Numerical integration result for the defect density as 
a function of r for 7 = 1 and 7 = 0.5. The parameter q is 
chosen such that the system crosses both the gapless phases 
in the two cases. For 7 = 1, we have used log-scale and slope 
of the straight line is 1/2 indicating a power- law decay of 
the defect density. For 7 = 0.5, we have used linear scale to 
accentuate the exponential fall of the defect density. 



diagram of the model for both 7=1 and 7 7^ 1 are 
shown in Fig. 1. For 7 7^ 1, excitation energy vanishes 
for the mode fc = tt at the phase boundary between an 
antiferromagnetic phase (AF) and a gapless phase (GPI) 
given by 2a = + (1 — 7)^. In the GPI phase, it 

is always possible to find a wave vector k for which 
vanishes. Similarly, vanishes for the mode /c = at 
the phase boundary given by 2a = -^/A^ (1 -f 7)^ which 
marks the boundary between the GPI and the second 
gapless phase (GPII). In GPII, both and vanish 
for some wave vector, so that there are four Fermi points 
in contrast to two Fermi points in GPI. The transition 
between GPI and GPII phases is a special quantum phase 
transition that involves doubling the number of Fermi 
points [23|. For 7 = 1, on the other hand, we arrive at a 
simplified form |A(fc)|2 = 2 + 2cosfc = 4cos2(fc/2). The 
phase boundaries between the antiferromagnetic phase 
and GPI phase, and the GPI and GPII phase are given by 
2a = A and 2a = VA^ -I- 4, respectively. It is noteworthy 
that for the case 7=1, the parameter |Ap vanishes 
at the AF-GPI phase boundary for k = tt and any a 
whereas in the anisotropic case, never does it vanish at 
the quantum transitions! 

In the quenching scheme we employ here, the scaled 
staggered magnetic field A is quenched as Xot/r from t — > 
--00 to +00 with a 7^ so that the system is swept across 
the quantum critical lines and the gapless phases. Let us 
set Ao = 1 for simplicity. As discussed already, so far as 
the dynamics is concerned, the term a cos k is irrelevant 
and the Schrodinger equation essentially reduces to a two 
level LZ problem. Note that in both the limits A ^• 
±00, the spins should be in a perfect antiferromagnetic 
orientation in the z-direction and wrongly oriented spins 
in the final state at t — > 00 (A — > 00) are the defects. 

Using the LZ transition formula we find that the 
probability of excitation in the final state is given by 
Pk = exp(-7r|A(fc)pT) = exp(-47rcos^(fc/2)T) for 7 = 1. 
In the adiabatic limit of r 00, only the modes close 
to fc ~ TT (for which |A(fc)p vanishes) contribute, so that 
Pk takes the simplified bompk = exp(— 7r(7r— k^r), and 



hence the defect density in the final state is given by 

I r 1 

n = — dk cxp{-n{n - k)'^T) ^ —=. (8) 
7^ Jo Vt 

For 7 ^ 1, on the other hand, A is minimum for k = tt 
on the AF-GPI phase boundary given by |Aop — |1 — 7^ 
so that the defect density 



n = — dfc cxp(— 7r|l -f- 7^ + 27COS fc|T) 
Jo 



We therefore, come across a special situation where the 
defect density decays exponentially with the rate r even 
though the system is swept across the critical lines and 
gapless phases as long as Ao 7^ 0, i.e., 7 7^ 1. We recover 
the power law scaling in Eq (8) for 7 = 1. As mentioned 
already, the scaling behavior shown in Eqs. (8) and (9) 
are same for all values of the scaled three spin interaction 
a. The defect density obtained by numerical integration 
of the Schrodinger equations (3) is shown in Fig 2. 

It is to be noted that the time evolution governed by 
Hk{t) is completely insensitive to the phase transitions 
generated by ac(fc). This fact can be understood also us- 
ing the following argument: the total number of Fermions 
for a mode fc is a constant of motion as far as the Landau- 
Zener dynamics is concerned. This is because the num- 
ber operator rik = a|,a/j -|- 6^,6^ commutes with Hk{t) 
for all k. In the initial (final) state A —00 ( +00), 
the expectation value of Uk is unity which as per the 
above argument is conserved throughout the dynamics. 
On the other hand, the three spin term gives rise to the 
gapless phases where both the energies and can 
become negative or both can become positive for some 
values of fc; the true ground state allows nk = 2 or 
in those cases. However, the instantaneous eigenstates 
of the time-dependent Hamiltonian will continue to have 
rifc = 1 as explained above and does not reach the true 
ground state. 
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Question remains whether the result presented here 
shouldpersist in the case of a general interacting system. 
In ref.Q, it was shown that the proof of the Kibble-Zurek 
scaling form does not require the system to be broken up 
into a product of two-level systems which can then be 
analyzed by Landau-Zener tunneling formula. The argu- 
ment only uses translational invariancc and some general 
scaling arguments, namely, the momentum dependence 
of the energies and the parameter dependence of the 
wave functions of those states. Extending the argument 
of ref.[H| in the present case demands that the Hamilto- 
nian should be decoupled into two parts; one time inde- 
pendent Hi{a) and the other time dependent H2{X{t)), 
where A = t/r as defined before, and Hi commutes with 
all the terms of H2{\). It can be shown that Hi{a) does 
not influence the dynamics except for a phase factor and 
the scaling form of the density of defects is then given 
by scaling form of the A-dependcnt part of the excitation 
energy (5a;2,fc(A) whereas vanishing of the total excitation 
energy SuJk = Stui^kict) + (5^2, fc (A) leads to the complex 
phase diagram with critical lines and gapless phases. If 
<5w2,fc(A) does not vanish, one expects an exponential de- 



cay of the defect density. On the other hand, if Suj2,k van- 
ishes for any mode k, then the density of defects decays 
as a powcrlaw with the quenching rate r. In the present 
work, we have used an example of an exactly solved sys- 
tem which satisfies the above conditions. Though pos- 
sible in principle , we believe that finding the example 
of such an interacting Hamiltonian which satisfies the 
above conditions and show an exponential decay of the 
defect density is a difficult and open problem. Our inter- 
est in this work is only to point out a special situation 
where one can find exponential decay of defect density 
even during passage through critical regions and we have 
illustrated the possibility using an exactly solvable spin 
chain with complicated interactions. 
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